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Abstract 

We  have  investigated  the  surface  chemistry  of  coadsorbed  hydrogen  and  halogen  atoms  on 
GeflOO),  produced  by  dissociative  chemisorption  of  HCI  and  HBr,  by  temperature-programmed 
desorption.  The  initial  sticking  probability,  Sq,  for  HCI  decreases  from  0.6  at  a  substrate 
temperature  of  270  K  to  0.05  at  400  K,  indicative  of  a  precursor  state  to  adsorption.  For  HBr  So  is 
constant  at  0.7  over  the  same  temperature  range.  A  fraction  /of  adsorbed  hydrogen  atoms  desorb 
associatively  as  H2  near  570  K,  while  the  remaining  (1-/)  H  atoms  recombine  with  adsorbed 
halogen  atoms  and  desorb  as  the  hydrogen  halide  (HX)  near  580  -  590  K.  The  activation  energies 
for  desorption  of  H2,  HCI,  and  HBr  are  all  approximately  40  kcal/mol.  For  both  HCI  and  HBr /  is 
0.7  at  low  initial  coverage  and  decreases  slightly  to  0.6  at  satiuation.  The  fraction /of  adsorbed 
halogen  atoms  left  on  the  surface  following  the  competitive  desorption  of  H2  and  HX  desorb  as  the 
dihalides  GeCl2  and  GeBr2  near  675  K  and  7 10  K,  respectively.  Desorption  of  H2,  HCI,  and  HBr 
occurs  with  near-first-order  kinetics,  similar  to  the  behavior  of  hydrogen  adsorbed  alone,  which  we 
attribute  to  preferential  pairing  induced  by  the  it  bond  on  unoccupied  Ge  dimers.  We  introduce  and 
solve  a  generalized  doubly-occupied  dimer  model  incorporating  competitive  pairing  of  H+H,  H-hX, 
and  X-)-X  on  Ge  dimers  to  explain  the  near- first  order  kinetics.  The  model  quantitatively  accounts 
for  both  the  desorption  kinetics  and  the  relative  yields  of  H2  and  HX  with  pairing  energies  of  *=  3 
kcal/mol.  Implications  of  the  present  results  for  surface  thermochemistry,  chemical  vapor 
deposition,  and  atomic  layer  epitaxy  of  Ge  and  Si  (100)2x1  surfaces  are  discussed. 
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I.  INTRODUCTION 


The  surface  chemistry  of  coadsorbed  hydrogen  and  halogen  atoms  plays  an  important  role  in 
the  growth  of  the  group  IV  semiconductors  silicon,  GejtSii-x  alloys,  and  diamond  by  chemical 
vapor  deposition  (CVD).  Surface  decomposition  of  chlorosilanes  on  silicon  surfaces  produces 
adsorbed  H  and  Cl,  and  in  fact  chlorosilane  CVD  has  historically  been  the  dominant  technology  for 
epitaxial  growth  of  silicon.*  Similar  techniques  have  recently  been  extended  to  growth  of  GexSii.jc 
alloys.2  H/Cl  surface  chemistry  has  also  been  implicated  in  determining  the  distribution  of 
crystallographic  orientations  that  remain  stable  during  silicon  growth  by  CVD.^  While  diamond 
film  growth  by  CVD  typically  involves  hydrogen  and  hydrocarbons,  several  authors  have  reported 
improved  results  using  alkyl  halide  precursors,*^-^  which  may  similarly  involve  coadsorbed 
hydrogen  and  halogen  atoms.  There  is  a  great  deal  of  current  interest  in  the  development  of 
methods  for  atomic  layer  epitaxy  (ALE)  techniques^  for  group  IV  materials,  and  hydrogen-halogen 
chemistry  offers  a  promising  approach.  Several  techniques  involving  alternating  cycles  of  a 
halogenated  precursor  and  hydrogen  have  been  proposed  for  ALE  of  silicon^*  *2  and  diamond.^ 
Besides  being  important  for  growth,  hydrogen-chlorine  chemistry  is  also  imponant  in  the  etching  of 
silicon.  *3 

In  addition  to  being  technologically  important,  coadsorbed  hydrogen  and  halogen  atoms  on 
semiconductor  surfaces  constitute  an  interesting  model  system  from  a  fundamental  point  of  view. 
Halogen  atoms  constitute  perhaps  the  simplest  adsorbates  on  semiconductor  surfaces  after 
hydrogen,  and  important  phenomena  such  as  ordering  and  intermixing  of  adsorbed  species  and  the 
branching  ratio  between  competing  reaction  channels  can  be  investigated  by  coadsorption  studies. 

To  date,  however,  only  a  handful  of  researchers  have  investigated  the  interaction  of  hydrogen 
halides*^20  or  dichlorosilane  (SiH2Cl2)^’^*  with  silicon  surfaces.  On  germanium,  only  very  early 
hydrogen-halogen  work  and  our  preliminary  results^^'^^  have  been  reported.  Further 

interest  in  hydrogen-halogen  surface  chemistry  stems  from  recent  observations  that  atomic 
hydrogen  can  readily  abstract  halogen  atoms  from  Si(100),25-27  which  is  directly  relevant  to 
ALE.5.9.11 

Much  more  attention  has  been  devoted  to  the  adsorption  of  hydrogen  or  halogen  atoms 
individually.  The  structures  of  adsorbed  hydrogen  on  Ge(  100)28-30  are  directly  analogous  to  those 
on  Si(l(X)),30-33  the  most  stable  species  being  the  monohydride,  with  one  hydrogen  atom  on  each 
dimerized  surface  atom.  Structures  analogous  to  the  monohydride  were  suggested  by  a 
photoelectron  spectroscopic  study  of  Cl  and  Br  on  Ge(  100)2x1,3'*  a  picture  which  is  supported  by 
the  corresponding  structures  of  C133.36  and  Br37  on  Ge(l  11)  and  of  Cl  on  Si(100)2xl.35^8-40 

The  mechanism,  kinetics,  and  dynamics  of  hydrogen  desorption  from  Si  and  Ge  (100)2x1 
surfaces  have  been  the  focus  of  a  great  deal  of  recent  effon  and  also  provided  much  of  the 
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motivation  for  the  present  study.  The  desorption  kinetics  are  nearly  first  order  in  the  hydrogen 
coverage  on  both  Si(100)4^‘^5  and  Ge(  lOO),''^  in  contrast  to  the  near-second-order  behavior  seen  on 
the  (111)  faces^2,47,48  and  on  metals.  The  desorption  kinetics  of  hydrogen  from  diamond  (100)  are 
also  approximately  first  order Hydrogen  molecules  desorbing  from  Si(l(X))  are  vibrationally 
hot50.5l  but  rotationally  cold;^®'^^  the  angular  distribution  is  peaked  toward  the  surface  normal^^ 
but  the  velocity  distribution  is  approximately  thermal.^^  The  available  experimental  evidence  seems 
to  favor  a  sequential^2,51  desorption  mechanism  between  H  atoms  paired  on  a  single 
dimer ."*2-46,49.5 1,54  xhe  near-first-order  kinetics  can  be  explained  by  preferential  pairing  of 
adsorbed  hydrogen, ^2^.54  which  is  driven  by  the  7t  bond^S  on  “unoccupied”  surface  dimers  and 
has  been  observed  directly  by  scanning  tunneling  microscopy  (STM).^"*  Fitting  the  deviation  from 
first-order  kinetics  at  low  coverage  to  the  predictions  of  our  doubly-occupied  dimer  model^2 
to  estimates  of  5-7  kcal/mol  and  =  5  kcal/mol  for  the  pairing  energies  of  hydrogen  on  Si(100)^345 
and  Ge(100),‘*^  which  is  also  a  measure  of  the  k  bond  strength  on  the  clean  surface.'*^’^  These 
estimates  of  the  pairing  energy  for  hydrogen  on  Si(l()0)  are  also  in  reasonable  agreement  with 
values  implied  by  recent  Si-H  bond  strength  calculations.32,33  However,  recent  high-level  quantum 
mechanical  calculations  of  the  desorption  dynamics  between  paired  hydrogen  atoms33. 56,57  on 
Si(  100)2x1  have  arrived  at  activation  energies  much  larger  than  the  57-58  kcal/mol  observed 
experimentally ."*2.44,45  The  apparent  discrepancy  between  theory  and  experiment  remains  to  be 
resolved. 

If  preferential  pairing  of  adsorbed  species  is  indeed  the  explanation  for  the  near-first-order 
desorption  kinetics  of  hydrogen  on  diamond,  Si,  and  Ge  (100)2x1  surfaces,  then  virtually  all 
recombinative  desorption  processes  between  mobile  adsorbates  on  group  IV  (100)2x1  should  be 
first  order43,54  since  the  driving  force  for  pairing,  the  dimer  tc  bond,  is  a  property  of  the  clean 
surface.  The  present  study  of  the  desorption  kinetics  of  HCl  and  HBr  from  Ge(lOO),  which  occurs 
in  competition  with  desorption  of  H2,  constitutes  the  first  test  cf  this  prediction  on  an  adsorbate 
other  than  hydrogen. 

More  limited  information  is  available  on  the  desorption  kinetics  of  adsorbed  halogens. 
Reactions  of  Cl2^*  and  Br2^9  with  Si(lOO)  have  been  investigated  by  temperature-programmed 
desorption  (TPD),  and  the  corresponding  reactions  with  Ge(lOO)  have  been  investigated  by  reactive 
scattering.^.6*  At  low  initial  halogen  coverages  or  high  reaction  temperatures  on  Si  or  Ge  (100) 
the  dihalide  is  the  dominant  desorption  (etch)  product,  while  tetrahalide  formation  becomes 
important  at  higher  initial  coverage  or  lower  reaction  temperatures. 

In  this  paper  we  report  the  adsorption  kinetics  of  HCl  and  HBr  on  Ge(lOO)  and  the  desorption 
kinetics  and  relative  yields  of  H2,  HX,  and  GeX2  (X=C1  or  Br),  as  determined  by  temperature- 
programmed  desorption  (TPD).  We  introduce  and  solve  a  generalized  doubly-occupied  dimer 
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model  that  incorporates  pairing  between  H,  X,  and  H+X  atoms  on  Ge-Ge  surface  dimers  and 
competitive  desorption  of  H2  and  HX,  and  show  that  the  new  model  quantitively  accounts  for  the 
data.  However,  the  inferred  pairing  energies  are  smaller  than  that  obtained  for  hydrogen  adsorbed 
alone.  The  results  are  discussed  within  the  context  of  estimates  of  the  thermochemistry  of  adsorbed 
H,  Cl,  and  Br  on  Ge  and  Si  surfaces. 

II.  EXPERIMENTAL 

Experiments  were  performed  in  an  ultrahigh  vacuum  (UHV)  chamber  that  is  described  in 
detail  elsewhere.^2,63  Briefly,  the  chamber  is  pumped  via  a  liquid-N2-trapped  diffusion  pump  and 
a  titanium  sublimation  pump,  and  is  equipped  with  LEED/ESDIAD  optics,  an  Auger  spectrometer 
(VSW  hemispherical  analyzer  HA- 100),  a  VG  SXP-400  quadrupole  mass  spectrometer  (QMS) 
with  a  water-cooled  shroud  mounted  on  a  translation  system,  a  calibrated  gas  doser,^  and  an  ion 
gun.  The  base  pressure  of  the  apparatus  was  1  x  10'^®  Torr  for  the  experiments  reported  here,  but 
some  of  the  experiments  were  performed  at  =  7  x  10*10  Torr  with  the  titanium  getter  saturated.  The 
sample  holder  was  mounted  on  a  rotatable  xyz  manipulator. 

The  Ge(lOO)  sample  was  cleaved  from  a  wafer  (Si-Tech,  Inc.)  cut  4-6°  off  the  (100)  plane 
towards  the  [011]  direction,  0.25-0.30  mm  thick,  n-type,  p  =  5-40  Q  cm,  into  a  rectangle  13.4  mm 
X  13.8  mm.  A  chromel-alumel  thermocouple  was  cemented  into  a  small  hole  (~  0.5  mm  dia.) 
drilled  near  one  edge  of  the  sample  using  Aremco  516  high-temperature  cement.  The  sample  was 
mounted  between  Ta-foil  clips  attached  to  a  Cu  block  and  could  be  heated  resistively  to  >  873  K  and 
cooled  with  liquid  nitrogen  down  to  153  K.  The  active  area  presented  to  the  doser  after  mounting 
was  1 1.5  mm  x  13.4  mm.  After  degreasing,  the  Ge  sample  was  placed  in  the  chamber  and  cleaned 
by  several  sputter- and-anneal  cycles  (jAr+  =  2-3  |iA  cm-^,  E\x*  =  5(X)  V,  Tanncai  =  850  K). 

Exposures  to  HCl  or  HBr  were  performed  by  rotating  the  sample  to  face  the  doser  and 
admitting  a  known  amount  of  gas  to  the  chamber  through  a  calibrated  aperture.  Computer 
calculations  of  the  gas  flux^^  show  that,  for  this  sample  geometry,  17  %  of  the  molecules  leaving 
the  doser  strike  the  sample,  yielding  a  flux  of  (4.1 1  ±  0.23)  x  10^^  P  S“^  where  M  is  the 
molecular  mass  (g  mol~0,  and  P  is  the  pressure  (Torr)  upstream  of  the  conductance-limiting 
orifice^  during  a  dose.  For  Ge(lOO),  one  monolayer  (ML)  =  6.23  x  10*^  atoms  cm**^,  and  is  used 
below  to  scale  both  doses  and  surface  coverages.  In  control  experiments,  the  sample  was  exposed 
to  atomic  hydrogen  by  backfilling  the  UHV  chamber  with  H2  to  pressures  of  2xl0'*-5xl0'^  Torr 
and  heating  a  coiled  W  filament  located  about  3  cm  from  the  sample  to  1700-1800  K. 

All  coverages  were  determined  by  TPD.  A  coverage  calibration  for  surface  hydrogen  was 
obtained  by  TPD,  detecting  H2  (m/e=2)  following  a  saturation  dose  of  H2S.  H2S  has  been  shown 
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to  adsorb  dissociatively  on  Ge(lOO)  as  H  +  up  to  0.5  ML  coverage  at  saturation,^^  and  yields 
exclusively  H2  and  GeS  as  desorption  products  upon  heating.62  We  used  the  same  coverage 
calibration  in  a  TPD  study  of  the  desorption  of  H2  from  Ge(lOO)  following  adsorption  of  atomic 
hydrogen,23.46  and  found  that  a  shoulder  at  525  K  in  the  TPD  peak  appeared  for  coverages  above 
0.96  ML.23  It  is  well  established  that  an  analogous  lower-temperature  TPD  peak  for  H2  desorbing 
from  Si(lOO)  appears  for  initial  coverages  greater  than  one  monolayer,^* •'^5,51,68  corresponding  to 
decomposition  of  dihydride  (SiH2)  species,  and  that  dihydride  species  on  Ge(lOO)  similarly 
decompose  at  a  lower  temperature  than  the  GeH  monohydride  state.^^  The  fact  that  the  onset  for 
the  appearance  of  the  shoulder  occurs  at  one  monolayer  therefore  supports  the  H2S  coverage 
calibration  described  above.  HX  (X  =  Cl  or  Br)  coverages  were  determined  by  assuming  that  both 
molecules  similarly  reach  a  saturation  coverage  of  0.5  ML  of  H  -t-  X.  This  last  assumption  is 
examined  critically  in  the  discussion  section  below. 

After  dosing,  when  the  background  pressure  returned  to  the  background  level,  the  sample  was 
rotated  to  face  the  entrance  slit  of  the  water-cooled  QMS  shroud,  at  a  distance  of  ==  0.5  cm.  The 
temperature  was  ramped  at  a  rate  of  2  K  s”*,  conffolled  by  a  Eurotherm  temperature  controUer,  and 
QMS  signals,  multiplexed  for  1-3  masses,  were  recorded  by  an  AT-compatible  personal  computer. 
We  found  that  the  most  reproducible  results,  particularly  with  HCl,  were  obtained  by  letting  the 
titanium  sublimation  pump  become  saturated  and  operating  at  a  chamber  pressure  of  7-8  x  lO’*® 
Torr.  The  reason  for  the  improved  reproducibility  is  presumably  that  the  chamber  pumping  speed 
remained  more  nearly  constant  so  that  integrated  TPD  peaks  were  accurately  proportional  to  the 
initial  coverage  regardless  of  the  recent  dosing  history.  Most  of  the  H2  and  HX  TPD  results 
presented  below  were  obtained  with  the  Ti  getter  saturated.  The  sublimator  was  used  when  GeX2 
desorption  data  was  collected  separately,  since  the  GeX2  species  are  expected  to  condense  on  the 
chamber  walls  regardless  of  the  presence  or  absence  of  a  getter. 


III.  RESULTS 

Three  desorption  products  were  observed  following  exposure  to  HX  (X  =  Cl  or  Br):  HX,  H2, 
and  GeX2,  as  shown  in  Figs.  1  and  2  for  HCl  and  HBr,  respectively.  In  both  cases  H2  desorption 
occurs  near  565  K  after  a  saturation  dose,  as  for  desorption  from  hydrogen  adsorbed  alone^^’"^  or 
following  decomposition  of  H2O  or  H2S.62  hCI  or  HBr  desorption  occurs  concurrently  with  H2 
desorption,  exhibiting  peak  temperatures  near  575  and  585  K,  respectively,  at  saturation  initial 
coverage.  Using  the  H2S  coverage  calibration  described  above  for  the  hydrogen  coverage  0.30 
ML  and  0.31  ML  of  hydrogen  atoms  desorb  as  H2  following  saturation  exposures  to  HCl  and  HBr, 
respectively.  Assuming  that  the  saturation  coverage  of  both  molecules  is  0.5  ML,  this  implies  that 
0.20  and  0. 19  ML  of  the  hydrogen  atoms  desorb  as  HX  for  HCl  and  HBr,  respectively.  The 
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dihalide  etch  products  GeCIi  and  GeBr2  desorb  at  higher  temperatures,  675  and  710  K, 
respectively,  following  saturation  HX  doses. 

The  magnitude  of  the  peak  desorption  temperature  for  HCl  and  HBr  provides  a  clear 
indication  that  the  adsorption  is  dissociative,  since  molecular  HCl  and  HBr  could  only  physisorb  to 
the  surface  and  would  be  expected  to  desorb  well  below  300  K.  However,  to  prove  that  the 
adsorption  is  dissociative  we  performed  control  experiments.  Adsorbed  Cl  or  Br  were  prepared  on 
the  surface  by  dosing  with  HCl  or  HBr,  respectively,  and  flashing  the  sample  to  approximately  635 
K.  This  procedure  removes  surface  hydrogen  as  H2  and  HX  (Figs.  1,  2),  leaving  adsorbed  halogen 
atoms,  and  avoids  the  complexity  of  separate  experiments  with  CI2  and  Br2.  The  sample  was  then 
rotated  to  face  the  W  filament  at  17(X)- 1 8(X)  K  and  given  a  nominal  dose  of  10  L  of  molecular  plus 
atomic  hydrogen.  The  sample  was  then  rotated  to  face  the  QMS  and  the  temperature  was  ramped. 
In  both  cases  an  HX  desorption  peak  was  observed  near  575  K.  Since  the  latter  experiment 
involved  reaction  between  absorbed  halogen  and  hydrogen  atoms  and  yielded  the  same  HX  TPD 
peaks  as  shown  in  Fig.  1  and  2  following  an  exposure  to  HX,  we  conclude  that  adsorption  of  HX  is 
indeed  dissociative. 

Initial  sticking  probabilities.  So,  for  HCl  and  HBr  were  obtained  from  the  slopes  of  coverage 
versus  exposure  data.  For  intermediate  HX  coverages,  the  coverage  of  hydrogen  desorbing  as  H2 
was  determined  as  the  ratio  of  the  integrated  TPD  peak  intensity  to  that  of  H2  following  a  saturation 
dose  of  H2S,  as  described  above.  The  coverage  of  hydrogen  desorbing  as  HX  was  determined  as 
the  ratio  of  the  integr  :ed  HX  peak  to  that  of  HX  following  a  saturation  dose  of  HX  multiplied  by 
the  inferred  coverage  of  hydrogen  desorbing  as  HX  at  saturation  (0.20  or  0. 19  ML  for  HCl  or  HBr, 
respectively).  The  sum  of  these  two  coverages  is  then  the  total  amount  of  hydrogen  initially 
chemisorbed  and,  by  stoichiometry,  the  initial  coverage  cf  dissociated  HX.  We  ran  a  number  of 
spot  checks  to  insure  that  the  amount  of  desorbing  H2  and  HX  following  a  null  dose  was  negligible. 
The  initial  sticking  probability  was  then  obtained  as  the  initial  slope  of  a  plot  of  total  HX  coverage 
versus  HX  exposure,  as  determined  using  the  calibrated  doser  described  in  the  preceding  section. 
The  results  for  Sq  as  a  function  of  substrate  temperature  are  summarized  in  Fig.  3.  Sq  for  HCl 
decreased  from  0.6  at  270  K  to  0.05  at  400  K,  while  that  of  HBr  was  constant  at  0.7  over  the  same 
temperature  range. 

Stoichiometry  implies  that  for  every  H2  molecule  that  desorbs,  the  two  X  atoms  remaining 
from  dissociative  adsorption  of  HX  must  desorb  as  GeX2.  The  desorption  peaks  for  H2  and  HX 
have  a  strong  overlap  while  those  of  GeX2  occur  at  higher  temperature.  Therefore,  the  branching 
ratio  (H2  +  GeX2  versus  HX  desorption)  may  be  usefully  described  by  the  ffaction/of  adsorbed 
hydrogen  atoms  which  desorb  as  H2-  The  observed  saturation  coverages  of  hydrogen  desorbing  as 
H2  of  0.30  and  0.31  ML  for  HCl  and  HBr,  respectively,  imply  that /  =  0.60  or  0.62  for  HCl  and 


HBr  from  the  assumed  saturation  HX  coverage  of  0.5  ML.  At  lower  initial  coverage  of  HCl  or 
HBr/ increases  slightly  to  0.66  and  0.68,  respectively,  in  the  low  initial-coverage  limit.  However, 
we  estimate  the  uncenainty  in  /  as  approximately  0.05,  so  the  increase  at  low  coverages  may  not  be 
real. 

TPD  spectra  for  H2  and  HX  as  a  function  of  the  initial  HX  coverage  are  shown  in  Figs.  4  and 
5  for  HCl  and  HBr,  respectively.  In  both  cases  H2  desorption  occurs  near  562  K  after  a  saturation 
dose,  and  the  peak  temperature  Tp  increases  to  about  580  K  at  initial  coverages  of  0.05  ML.  For 
HCl  and  HBr  desorption  Tp  increases  from  578  and  583  K,  respectively,  at  saturation  initial 
coverage  to  590  and  595  K  at  initial  coverages  of  =0.05  ML.  The  weak  dependence  of  Tp  on  initial 
coverage  for  both  H2  and  HX  desorption,  together  with  the  asymmetric  peak  shapes  (Figs.  1,2,4, 
and  5),  indicates  near-first-order  kinetics.  The  dependence  of  Tp  for  both  H2  and  HX  on  the  initial 
HX  coverage  is  summarized  in  Figs.  6  and  7  for  HCl  and  HBr,  respectively,  together  with  fits 
made  using  the  model  described  in  the  next  section. 

We  have  not  investigated  the  coverage  dependence  of  the  GeCl2  and  GeBr2  desorption 
products  in  detail  because  the  signal-to-noise  ratio  in  the  data  was  significantly  worse  than  that  for 
H2  and  HX  and  because  the  modification  of  the  surface  structure  implied  by  the  Ge  atom  renxjval 
associated  with  desorption  complicates  the  interpretation  of  the  results.  We  found  that  Tp  for  GeCl2 
was  constant  at  675  K  to  within  10  K  for  initial  HCl  coverages  between  0.03  and  0.5  ML.  For 
GeBr2  Tp  increased  from  7 10  K  at  saturation  coverage  to  765  K  at  an  initial  HBr  coverage  of  0.05 
ML.  The  desorption  behavior  of  GeCl2  is  indicative  of  first-order  kinetics  while  that  of  GeBr2 
suggests  a  kinetic  order  intermediate  between  1  and  2. 

IV.  COMPETITIVE  PAIRING  MODEL 

The  near-first-order  desorption  kinetics  of  HX  implied  by  the  weak  dependence  of  on  initial 
coverage  are  qualitatively  similar  to  the  desorption  behavior  of  hydrogen  on  Ge(l(X))2xl'*6  and  on 
Si(l(X))2xl.'^l‘^5  ^5  previously,  for  hydrogen  adsorbed  alone,^^-'^  we  attribute  the  near-first-order 
kinetics  to  preferential  pairing  of  adsorbed  atoms  on  surface  dimers  driven  by  the  weak  k  bonds  on 
clean  surface  dimers.  The  presence  of  two  types  of  atoms  (H,  X),  three  types  of  pairing  (H+H, 
H+X,  X-t-X),  and  three  desorption  channels  necessitates  a  generalization  of  the  doubly-occupied 
dimer  model  that  we  developed  to  explain  the  desorption  kinetics  of  hydrogen  from  Si(l(X))2xl'*3 
and  is  described  below. 

Consider  a  surface  with  M  (2x1)  dimer  sites,  each  of  which  can  be  empty,  or  singly-occupied 
by  H  or  X,  or  doubly-occupied  by  H-t-H,  H+X,  or  X+X.  Suppose  that  2Nh  hydrogen  atoms  and 
2Nx  halogen  atoms  are  chemisorbed  on  the  surface  and  that  Nn+Nx  <  M,  so  that  the  (2x1) 
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reconstruction  is  undisturbed.  For  the  present  we  assume  that  the  dinars  are  non-interacting  as  in 
the  original  model.^^  We  consider  the  effects  of  interactions  between  dimers  for  H  on  Si(  100)2x1 
elsewhere.^^'^0  por  simplicity  let  us  assume  that  in  the  (degenerate)  ground  state  of  the  system,  all 
2A^H  hydrogen  atoms  are  paired  on  jVh  dimer  sites  and  all  2Nx  halogen  atoms  are  paired  on  Nx 
dimer  sites.  Identical  results  are  obtained  if  the  ground  state  is  instead  assumed  to  consist  of  paired 
H+X  atoms  plus  paired  H  or  X  atoms,  whichever  species  is  present  in  excess.  Let  Chh  and  exx  be 
the  energies  required  to  unpair  the  hydrogen  or  halogen  atoms,  respectively,  that  is,  to  move  a  H  or 
X  atom  from  a  doubly-occupied  to  an  empty  dimer,  creating  two  singly-occupied  dimers,  as 
indicated  in  Eqs.  (1)  and  (2); 


H  H 

Ge-<je 

^  \ 

X  X 

be-cfe 


Ge=Ge 


Ge=Ge 


H 

\  • 
Ge-Ge 


^e-Ge 


Suppose  there  are  «h  and  nx  such  unpairings,  respectively,  in  the  system.  The  H  and  X  atoms  on 
singly-occupied  dimers  produced  by  unpairings  of  HH  and  XX  dimers  can  pair  up  to  form  an  HX 
doubly-occupied  dimer,  suppose  there  are  whx  such  re-pairings.  The  energy  gained  by  re-pairing  of 
H-t-X  is  CHx.  as  indicated  in  Eq.  (3): 


H 

^e-Ge 

^  \ 


^e-Ge 

7  ^ 


H  X 

^e-cfe 

f  \ 


Define  mHH.  'whx.  and  mxx  to  be  the  numbers  of  doubly-occupied  dimers  with  paired  H+H,  H+X, 
or  X+X  atoms,  respectively,  mn*  and  mx*  to  be  the  numbers  of  singly-occupied  dimers  with  H  or 
X  atoms,  and  m**  to  be  the  number  of  unoccupied  dimers.  It  is  readily  seen  that 

"iHH  =  Nil-  nu 

tnxx  =  A'x  -  nx 

"JHX  =  nHX  ... 

(4) 

mu*  =  2nn  -  rtHX 
mx*  =  2nx-nHX 
m**  =  M  -  Nw  -  Nx  -  nyi  -  nx  +  njfx 

The  partition  function  Z  for  this  non-interacting  lattice  gas  model  is  only  marginally  more 
complicated  than  that  for  a  one-component  system:^^ 


nH.nX.^HX 


jiff  +  "IH*  ^^nx*)  J  13  1  /C\ 

T  '  - j - : - j - T  exp[  -J3  (nwQ/w  +  nx^X  *  n«X^x)J  (5) 

H  H  •  niH  X  •  ^nx  x  •  nzH  * !  mx  • '  nt*  • ! 


where  s  1/IcqT.  The  combinatorial  factor  arises  from  the  number  of  ways  of  arranging  the  six 
type  of  dimers,  and  the  2^"^* + /r*i*  +  mx*)  factor  takes  account  of  the  fact  that  H+X-occupied  dinars 
and  both  H*  and  X*  singly-occu^i^d  dimers  can  each  be  arranged  in  two  different  ways. 
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The  equilibrium  pop’  udons  of  the  various  dimer  arrangements  are  obtained  by  maximizing 
the  partition  function.  Approximating  InZ  by  the  logarithm  of  the  maximum  summand,  using  Eqs, 
(4)  to  replace  the  numbers  of  dimers  by  the  numbers  of  excitations  «h.  'tic.  and  /ihx.  and  setting  the 
partial  derivatives  of  InZ  derivative  with  respect  to  «h.  nx.  and  /jhx  equal  to  zero,  we  obtain  three 
equations  in  the  three  unknowns: 


4exp(-/3eHH)  = 
4exp(-/3exx)  = 
Icxpi-peux)  = 


_ (2/iH  -  hhxF _ 

(A/h  -  -  A'h  -  i'Vx  -  «h  -  «x  +  «Hx) 

_ (2nx  -  _ 

{N\  -  nx)(M  -  A/h  -  A/x  -  -  nx  +  «Hx) 

(2/»h  -  nHx)(2/tx  -  hhx) 

/2HX  (iW  -  A/h  -  A/x  -  riH  -  «X  +  WHx) 


(6a) 

(6b) 

(6c) 


We  now  introduce  the  coverages  of  hydrogen  and  halogen  atoms,  expressed  with  respect  to 
the  number  of  dangling-bond  sites  and  coverages  of  the  various  types  of  doubly-occupied  and 
singly-occupied  dimers: 


= 

NhIM 

&K  = 

NxIM 

= 

"Jhh/A/ 

&XX  = 

mx\IM 

0m  = 

mwx/Af 

= 

m\{*IM 

&x*  = 

mx*IM 

&*•  = 

m**IM 

Re-introducing  the  dimer  numbers  defined  in  Eq.  (4)  into  Eqs.  (6)  and  introducing  the  coverages 
defined  in  Eq.  (7),  we  obtain  three  rather  simple  equations: 


4exp(-/3eHH) 

_  0H* 

(8a) 

6>hh^* 

4  c\p(-j5exx) 

II 

* 

(8b) 

&XX0*- 

2  exp(-/JeHx) 

_  &H*0X* 

(8c) 

<9hx^* 


0 


Eqs.  (8)  may  be  solved  most  conveniently  by  choosing  as  the  three  independent  variables  6^, 
6^]*,  and  It  is  easily  seen  that  the  remaining  dimer  coverages  can  be  expressed  in  terms  of 
these  three  together  with  the  total  hydrogen  and  halogen  coverages,  ©h  and  respectively,  by 


+  (9a) 

^tx  =  ^  +  ^tx)  (9b) 

O**  =  I  -  6^1 '  ©x  -  ^*)  (9c) 

Eq.  (8c)  may  then  be  trivially  solved  to  obtain  0hx  in  terms  of  Oh*  and  0^*: 

Oqx  =  1-  exp(/3£Hx)  (10) 

2  O*. 


where  O**  is  given  by  Eq.  (9c).  Unfortunately,  Eqs.  (8a)  and  (8b)  do  not  simplify  further,  leaving 
two  coupled  equations  in  the  unknowns  Oh*  and  Ox*  which  must  be  solved  numerically. 

Solution  of  Eqs.  (8-10)  yields  the  coverages  of  each  possible  combination  of  species  adsorbed 
on  dimers  for  a  specified  total  hydrogen  and  halogen  coverage  at  a  given  temperature.  All  that 
remains  to  be  able  to  simulate  the  TPD  data  is  to  postulate  a  relationship  between  the  dimer 
coverages  and  the  desorption  rates. 

We  postulate,  as  before, that  desorption  occurs  between  atoms  paired  on  a  single  dimer, 
which  in  this  case  includes  H+X  as  well  as  H+H.  The  desorption  rates  of  H2  and  HX  are  then 
given  by  and  A:Hx0kx.  rrspectively,  where  ^hh  and  ^hx  are  the  rate  constants  for 

desorption  of  H2  and  HX.  respectively.  We  neglect  desorption  of  GeX2  since  it  occurs  at  higher 
temperature  and  we  seek  to  model  only  the  competitive  desorption  of  H2  and  HX.  The  kinetic 
equations  for  the  surface  coverage  during  a  TPD  experiment  are  then 

=  -  -  ^Hx^kx  (11a) 

=  -^HX0kx  (11b) 

We  assume  that  the  preexponential  factor  for  ^hh  is  2x10*^  sec'^  as  was  found  for  H2  desorption 
from  81(100)  by  Hofer  et  who  took  H-atom  pairing  into  account  in  the  analysis,  and  a  value  of 
lOl^  sec*^  for  the  preexponential  factor  for  knx-  As  discussed  elsewhere,^  the  dependence  of  TPD 
peak  temperature  on  coverage  is  quite  insensitive  to  the  value  of  the  preexponential  factor.  The  TPD 
peak  temperatures  for  H2  and  HX  at  high  initial  coverage  then  imply  activation  energies  of 
approximately  42  and  38  kcal/mol  for  /:hh  and  A:hx.  respectively. 


dt 

dt 
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Simulated  TPD  experiments  at  initial  coverages  of  ^  ^  between  0  and  O.S  ML  began  at  a 

surface  temperature  of  450  K,  which  was  then  ramped  linearly  at  2  K  s' *  as  in  the  experiments.  At 
each  temperature,  the  equilibrium  surface  dimer  coverages  (Eqs.  7)  were  calculated  by  solving  Eqs. 
(8-10)  numerically,  evaluating  the  rate  constants  /:hh  and  ^hx.  and  calculating  the  time  derivatives  of 
^  and  (Eqs-  1  !)•  The  use  of  equilibrium  dimer  coverages  assumes,  of  course,  that  surface 
diffusion  is  sufficiently  fast  relative  to  desorption  for  quasiequilibrium  to  be  maintained  even  as 
paired  hydrogen  and  H+X  are  removed  by  desorption  (Eq.  1 1).  The  dynamic  surface  coverages 

and  Gxit)  were  generated  in  tandem  using  a  4th  order  Rungc-Xutta  integration  scheme.  The 
desorption  rates  for  H2  and  HX  were  stored  at  each  time  step.  The  time  increment  was  chosen  so 
that  the  surface  temperature  increases  by  1  K  at  each  step,  which  we  found  to  yield  converged 
results.  \s  H2  desorbs  the  surface  becomes  increasingly  enriched  in  X  relative  to  H  (6^  <  Sx:)- 
After  each  run,  the  peak  temperatures  Tp  for  H2  and  HX  and  the  fraction  /  of  hydrogen  desorbing 
as  H2  were  evaluated.  The  pairing  energies  were  then  determined  by  performing  a  least-squares  fit 
between  measured  and  predicted  values  of  the  peak  temperature  and  H2  yields.  The  activation 
energies  for  ^hh  and  ^hx.  £^a.HH  and  f'a.HX.  respectively,  were  also  allowed  to  vary  so  as  to  allow 
for  small  shifts  in  Tp,  with  the  rationale  that  coadsorbed  halogen  atoms  may  slightly  perturb  the 
desorption  kinetics. 

Excellent  agreement  between  experiment  and  model  predictions  for  both  HCl  and  HBr  was 
obtained  by  assuming  equal  values  for  Chh.  ^hx*  and  €xx-  For  HCl  best  agreement  was  obtained 
from  £Hh  =  =  £xx  =  2.7  kcal/mol,  £a,HH  =  ^1-4  kcal/mol,  and  £a,HX  =  38.1  kcal/mol.  With 

these  parameters  the  fraction  /  of  surface  hydrogen  desorbing  as  H2  is  68%  at  0.5  ML  initial 
coverage  and  increases  very  slightly  at  lower  coverage  (69%  at  0.05  ML).  The  model  predictions 
for  Tp  are  plotted  together  with  the  experimental  results  in  Fig.  6.  The  corresponding  best 
parameter  values  for  HBr  are  £hh  =  ^HX  =  ^xx  =  2.8  kcal/mol,  £a,HH  =  41.7  kcal/mol,  and  £a,HX  = 
38.5  kcal/mol.  The  predicted  values  of  7p  for  H2  and  HBr  are  shown  in  Fig.  7  with  the  data.  For 
HBr /  is  calculated  as  69%  at  0.5  ML  initial  coverage  and  again  increases  very  slightly  at  lower 
coverage  (70%  at  0.05  ML).  Agreement  between  experimental  and  model  results,  for  Tp,/,  and 
their  coverage  dependences,  is  quite  good  for  both  HCl  and  HBr.  The  results  of  fits  obtained  using 
different  choices  of  parameters  are  discussed  below. 


V.  DISCUSSION 
A.  Adsorption 

The  adsorption  kinetics  for  HCl  and  HBr  on  Ge(lOO)  are  consistent  with  the  paradigm  we 
have  proposed^3 ,62,71  for  surface  chemistry  on  group  IV  (l(X))2xl  surfaces;  the  dimerized  surface 
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atoms  are  analogous  to  highly  strained  digermenes  (compounds  containing  Ge-Ge  double  bonds), 
disilenes  (Si=Si  compounds),  or  olefins  on  GeflOO),  Si(lOO)  or  diamond  (100),  respectively.  If 
reaction  occurs  on  a  single  surface  dimer,  adsorption  and  desorption  are  analogous  to  molecular 
addition  and  elimination  reactions,  respectively.  Disilene  and  digermene  compounds  have  been 
known  for  more  than  a  decade,  being  synthesized  for  the  first  time  in  1981^2  and  1982,’^3 
respectively.  The  it  bonds  in  these  molecules  are  considerably  weaker  than  in  olefins  and  readily 
undergo  addition  reactions  with  HCl,  halogens,  water,  and  alcohols.  Several  excellent  reviews  of 
disilene  and  digermene  chemistry  are  available. 

The  high  sticking  probabilities  for  HCl  and  HBr  on  Ge(lOO)  (Fig.  3)  imply  negligible 
activation  energies  for  adsorption  of  both  molecules.  The  absence  of  an  activation  barrier  suggests 
that  adsorption  occurs  through  four-center  transition  states  on  a  single  surface  dimer,^2  jg  that  Ge- 
H  and  Ge-X  bonds  can  begin  to  form  before  the  H-X  bond  is  broken: 


H-X 

Ge=Ge. 


H-X 


1$ 


Gc^Ge  , 


Ge-Ge 

f  \ 


(12) 


Simple  geometric  considerations  suggest  that  dissociation  of  HX  across  a  single  dimer  occurs  more 
readily  than  between  dangling  bonds  on  adjacent  dimers.  The  Ge=Ge  dimer  bond  length  on  the 
clean  surface  may  be  estimated  as  2.41  A  from  the  measured  value  of  its  parallel  component  (2.34 
A)^5  and  the  buckling  angle  (14°)  calculated  by  high-level  density  functional  theory.^^  The  bond 
lengths  in  HCl  (1.27  A)  and  HBr  (1.41  A)  are  shorter  than  the  Ge=Ge  dimer  bond,  and  therefore 
bond  strain  will  be  present  in  the  transition  state  (Eq.  12).  Ge  atoms  in  adjacent  dimers  in  the  same 
row  are  separated  by  a  considerably  larger  distance  (3.99  A).  The  HX  bond  would  need  to  be 
nearly  broken  before  Ge-H  and  Ge-X  bonds  on  adjacent  dimers  could  be  formed,  which  should  be 
accompanied  by  a  substantial  activation  energy.  Dissociative  adsorption  on  a  single  dimer  therefore 
seems  more  likely.  The  Ge-Ge  dimer  bond  length  is  intermediate  between  the  bond  length  in 
molecular  digermenes,  2.21-2.35  A,^3.77  (he  single  bond  length  of  2.44  A  in  bulk  germanium, 
as  would  be  expected  for  a  highly  strained  double  bond.  Ge  surface  dimers  should  therefore  be 
even  more  reactive  than  molecular  digermenes.  Silicon  and  diamond  surface  dimers  likewise 
should  be  substantially  more  reactive  than  the  analogous  disilene  and  olefin  compounds.  Four- 
center  transition  states  for  molecular  addition  reactions  with  7t-bonded  silicon  and  germanium 
compounds  have  been  proposed  since  the  early  1960’s,^^*^^  and  therefore  it  is  reasonable  to 
suppose  that  similar  transition  state  configurations  occur  during  dissociative  adsorption  on  7t- 
bonded  surface  dimers.  Although  we  have  represented  the  transition  state  as  being  nearly 
symmetric  for  ease  of  visualization,  the  HX  moiety  is  likely  to  be  tilted  and  shifted  toward  one  end 
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of  the  dimer,  with  the  Ge-X  distance  being  longer  than  Ge-H  distance,  based  on  the  analogous 
transition  state  calculated  for  addition  of  HCl  to  H2Si=SiH2.^ 

The  decrease  of  So  for  HCl  with  surface  temperature  is  analogous  to  the  behavior  seen  for 
H2O  on  Ge(iOO)^2,80  and  suggests  the  existence  of  a  mobile  physisorbed  precursor  state.^^  The 
classical  precursor  model  of  Kisliuk,^!  which  assumes  thermalization  of  physisorbed  molecules  in 
a  “precursor”  state,  predicts  a  temperature-dependent  initial  sticking  probability  given  by 

SoiT)  =  - CE - ^ -  (13) 

1  +  Vd/Vc  exp[-(£;i  -  Ec)/RT] 

In  Eq.  (13),  the  trapping  probability  a  is  the  probability  that  the  incident  molecule  becomes 
“trapped”  into  the  physisorbed  precursor  state,  Vd  andVg  are  the  preexponential  factors  for 
desorption  and  chemisorption,  respectively,  from  the  precursor  state,  and  £d  and  Eq  arc  the 
corresponding  activation  energies.  In  the  absence  of  a  direct  measurement  of  a  for  HCl  on 
Ge(lOO),  we  simply  assume  a  =  1 .  The  best  fit  to  the  data  is  then  obtained  with  values  of  7  x  10^ 
for  Vd/Vc  and  6.3  kcal/mol  for  Ed-^^c-  The  fit  is  shown  together  with  the  data  in  Fig.  3.  The 
modest  value  of  Ed-^^c  is  typical  for  precursor  adsorption  systems.  Only  a  limited  significance 
should  be  ascribed  to  these  precursor  model  parameters  due  to  the  uncertainty  in  the  data  and  the 
oversimplification  of  the  adsorption  dynamics  by  the  model. 

While  the  temperature  independence  of  So  for  HBr  could  be  explained  by  a  precursor  model 
with  Ed  fortuitously  equal  to  Ec,  we  suggest  that  a  more  likely  explanation  is  that  Ec,  the  activation 
energy  for  dissociative  chemisorption  for  physisorbed  HBr,  is  approximately  equal  to  zero.  In  this 
case  incoming  molecules  would  never  thermalize  into  a  physisorbed  state  but  would  simply 
dissociate  at  an  adsorption  site  once  enough  translational  energy  had  been  dissipated  to  the  lattice  for 
the  molecules  to  become  localized.  Sq  would  simply  equal  a,  where  now  a  is  interpreted  as  the 
fraction  of  incoming  molecules  that  lose  enough  translational  energy  to  the  surface  during  the  initial 
collision  to  undergo  subsequent  “bounces.” 

Why  is  So  for  HBr  on  GeflOO)  larger  and  Ec  apparently  smaller  than  the  corresponding 
quantities  for  HCl?  The  simplest  explanation  would  be  a  larger  heat  of  adsorption  for  HBr,  since 
the  more  exothermic  of  a  pair  of  similar  chemical  reactions  typically  has  the  smaller  activation 
energy .82  However,  estimates  based  either  on  bond  strengths  in  analogous  molecules  or  on  the 
activation  energies  for  desorption  of  HCl  and  HBr  lead  to  the  conclusion  that  the  heats  of  adsorption 
of  HCl  and  HBr  on  GeflOO)  are  very  nearly  equal,  as  discussed  below.  In  addition,  a  direct 
calorimetric  measurement  of  the  heat  of  adsorption  on  polycrystalline  Ge  powder  yielded  nearly 
identical  values  for  HCl  and  HBr.  ^<5  suggest  an  alternate  simple  explanation  for  the  difference 
in  5o,  motivated  by  the  hypothetical  transition  state  (Eq.  12).  Because  the  H-Br  bond  is  longer  than 
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the  H-Cl  bond,  less  distonion  is  necessary  to  achieve  the  significant  orbital  overlap  with  the  orbitals 
on  the  Ge  dimer  atoms  necessary  for  a  low-energy  transition  state.  In  addition,  the  force  constant 
for  the  H-Br  bond  is  20%  smaller  than  that  for  HCl,^^  ^  that  the  necessary  distortion  costs  less 
energy.  Further  experiments  on  the  reactivity  of  homologous  molecules  with  Ge(  100)  and  Si(lOO) 
are  needed  to  determine  whether  this  line  of  reasoning  is  capable  of  predicting  trends  in  sticking 
probabilities  correctly. 

Independent  evidence  for  dissociative  adsorption  across  single  dimers  on  Si(  100)2x1, 
consistent  with  the  adsorption  model  depicted  in  Eq.  (12),  comes  from  recent  investigations  of 
water  adsorption  by  scanning  tunneling  microscopy  (STM).*^'85  Two  recent  smdies  have  provided 
clear  evidence  that  many  or  most  of  the  defects  previously  observed  on  Si(lOO)®^  and  presumed  to 
be  intrinsic  are  in  fact  products  of  water  adsorption.^'*^  jhe  most  prevalent  of  these  “dark”  sites 
are  dimers,  suggesting  that  dissociative  adsorption  as  H+OH  on  a  single  dimer  is  the  most 
common  adsorption  mechanism.  The  STM  images  also  provide  evidence  for  adsorption  on 
adjacent  dimers,^'^^  including  isolated  dangling  bonds  of  low  reactivity  that  were  left  over  when 
most  pairs  of  adjacent  dangling  bonds  became  saturated.^  We  note  that  dissociation  fragments 
produced  on  a  single  dimer  could  become  localized  on  adjacent  dimers  even  at  low  temperature  if 
exothermic  bond  formation  during  the  chemisorption  process  causes  a  transient  local  heating, 
enabling  one  species  to  hop  to  an  adjacent  dimer. 

We  believe  that  the  applicability  of  the  molecular  paradigm  (Eq.  12)  for  the  surface  reactivity 
of  Ge(lOO)  and  Si(lOO)  should  be  examined  in  further  detail  to  determine  whether  it  can 
consistently  rationalize  reactivity  trends.  If  it  can,  predictive  capability  should  be  possible,  which 
could  be  useful  to  a  wide  variety  of  semiconductor  growth  and  processing  applications. 

B.  Desorption  of  H2,  HCI,  and  HBr  and  thermochemical  estimates 

By  detailed  balance,  if  adsorption  of  HCI  and  HBr  occurs  on  a  single  Ge-Ge  surface  dimer 
(Eq.  12)  then  the  reverse  process,  desorption,  will  occur  between  paired  H+Cl  or  H+Br  on  a  single 
dimer,  as  was  explicitly  assumed  in  the  competitive  pairing  model  described  in  section  FV.  A  local 
picture  for  the  desorption  process  is  strongly  supported  by  the  insensitivity  of  the  kinetics  of  H2 
desorption  from  Ge(l(X))  to  coverages  of  Cl,  Br,  O,  or  S  as  high  as  0.5  ML  (produced  by 
dissociative  adsorption  of  HCI,  HBr,  H2O,  or  H2S)^^’^2  This  insensitivity  to  coadsorbed  species 
contrasts  strongly  with  desorption  kinetics  on  metal  surfaces,  where  bonding  is  much  more 
delocalized,  which  are  often  strongly  perturbed  by  even  a  few  percent  of  a  monolayer  of  coadsorbed 
species. 
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The  central  conclusion  of  this  paper  is  that  preferential  pairing  of  adsorbed  H  and  X  must  be 
postulated  in  order  to  account  for  the  near- first-order  desorption  kinetics.  We  have  shown  that  the 
competitive  pairing  model  accurately  accounts  for  both  the  coverage  dependence  of  Tp  and  the 
relative  yields  of  H2  and  HX  with  pairing  energies  of  H+H,  H+X,  and  X-i-X  of  about  3  kcal/mol. 

To  assess  the  robustness  of  the  overall  conclusion  and  the  pairing  energies  we  have  performed  a 
number  of  additional  calculations.  The  sensitivity  of  the  fits  of  the  HCl  data  to  the  assumed  pairing 
energies  is  illustrated  in  Fig.  8.  Here  the  common  pairing  energies  were  fixed  at  1, 2.7,  and  5 
kcal/mol  and  the  activation  energies  for  desorption,  £a.HH  and  £a,HX.  were  adjusted  slightly  for  the 
best  fit  to  the  data.  It  is  clear  that  significantly  larger  or  smaller  pairing  energies  yield  a  much  poorer 
fit  to  the  data.  We  estimate  an  uncertainty  of  ±1  kcal/mol  in  the  common  effective  pairing  energy. 

Although  our  data  are  insufficient  to  establish  how  nearly  equal  the  pairing  energies  £)hh. 
and  £xx  are  to  one  another,  we  can  say  that  £hx  must  be  comparable  in  magnitude  to  eHH-  In  other 
words,  the  data  cannot  be  fit  by  assuming  preferential  pairing  of  hydrogen  only.  Shown  in  Fig.  9 
are  the  results  of  fits  to  the  HCl  data  made  assuming  that  one  or  both  of  £Hx  and  exx  are  zero.  If 
both  £HX  and  exx  are  set  equal  to  zero,  a  noticeably  worse  fit  to  the  dependence  of  Tp  on  6b  is 
obtained  (dashed  curves).  However,  the  most  serious  deficiency  of  a  fit  omitting  preferential 
pairing  of  HX  (i.e.,  Chx  =  0)  is  that  the  coverage  independence  of  Tp  for  HX  can  only  be 
reproduced  if  nearly  all  the  hydrogen  desorbs  as  H2  (i.e.,/*  I)  so  that  HX  desorption  becomes 
pseudo  first  order.  An  improved  fit  to  the  Tp  data  incorporating  preferential  pairing  only  of 
hydrogen  can  only  be  obtained  with  even  worse  agreement  with  the  H2  yields.  We  thus  conclude 
that  Chx  is  comparable  in  magnitude  to  £hh-  A  less  definitive  statement  is  possible  whether  £xx  is 
comparable  in  magnitude  to  chh  and  f^x-  The  fit  to  the  Tp  data  assuming  £Xx  =  0  and  equal  values 
of  and  £!«  (solid  curves)  is  better  but  still  inferior  to  the  case  where  all  the  pairing  energies  are 
equal  (Fig.  6).  The  predicted  yields  in  this  case  are  good,  as  is  the  dependence  of  Tp  on  6b  for  HCl. 
However,  the  predicted  values  of  Tp(6^)  for  H2  exhibit  a  much  stronger  curvature  than  that  seen 
expe  mentally  or  that  predicted  by  the  competitive  pairing  model  with  ^  We  thus 
conclude  that  £Xx  is  also  nonzero,  but  cannot  make  a  definitive  statement  as  to  how  close  in 
magnitude  it  is  to  the  other  pairing  energies. 

It  is  not  clear  whether  any  inference  about  the  value  of  exx  can  be  made  from  the  desorption 
kinetics  of  GeX2.  Each  desorption  event  removes  one  Ge  atom  from  the  surface,  presumably 
disrupting  the  dimer  structure.  This  will  affect  the  driving  force  for  pairing  and  also  the  surface 
mobilities,  as  diffusion  along  the  dimer  rows  is  likely  to  be  more  facile  that  diffusion  perpendicular 
to  them.  Annealing  of  the  Ge  surface  will  also  occur  at  these  temperatures,  but  in  the  absence  of 
evidence  that  annealing  of  the  etch  damage  is  fast  in  comparison  to  desorption  we  believe  that  an 
analysis  of  the  desorption  kinetics  of  GeX2  based  on  an  ideal-dimer  model  is  not  justified. 
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Finally,  we  have  also  fit  the  data  using  a  phenomenological  desorption  model  wherein  pairing 
is  not  treated  explicitly,  as  described  in  more  detail  elsewhere.^^  Rather  then  assuming  that  the 
desorption  rates  of  H2  and  HX  are  and  respectively,  we  set  these  quantities  equal 

to  atid  The  rate  constants  km  and  i^HX  have  tlic  sanM  meaning  as  in  the 

competitive  pairing  model,  but  pairing  is  treated  indirectly  via  the  phenomenological  reaction  orders 
n  and  m.  With  these  substitutions,  the  kinetic  equations  for  the  surface  coverage  during  a  TPD 
experiment  become 

at 

^  (14b) 

The  simulation  of  TPD  experiments  was  otherwise  carried  as  for  the  competitive  pairing  model. 
The  preexponential  factor  for  HX,  Ahx.  was  allowed  to  vary  along  with  n  and  m,  while  the 
prexponential  factor  for  H2  was  fixed  as  2  x  10*5  s‘*  and  both  £a,HH  and  £a.HX  were  set  equal  to  42 
kcal/mol.  Fits  roughly  comparable  in  quality  to  those  of  the  competitive  pairing  model  were 
obtained  for  the  coverage  dependence  of  Tp  for  both  H2  and  HX  and  the  relative  desorption  yields. 
For  HCl  the  best-fit  parameters  were  «  =  1.2,  m  =  0.4,  and  Ahx  =  5  x  lO**^  s**.  For  HBr  the 
corresponding  values  were  n  =  1.35,  m  =  0.5,  and  Ahx  =  4  x  10*^  s'*. 

Although  the  underlying  physics  of  the  phenomenological  model  are  not  explicit,  the 
conclusion  from  the  fitting  is  the  same  as  with  the  competitive  pairing  model:  correlation  of  both 
H-fH  and  H+X  is  necessary  to  fit  the  data.  The  parameter  n  is  the  kinetic  order  for  H2  desorption. 
The  vales  of  n  close  to  1  reflect  the  nonrandom  distribution  of  H  atoms  and  concomitant  near-first- 
order  desorption  kinetics.  A  good  fit  could  not  be  achieved  with  m  =  I,  which  would  imply 
random  distributions  of  H  and  X:  0hx  =  Values  of  m  less  than  one  correspond  to 

correlation  of  adsorbed  H  and  X  atoms  (m  =  0  would  imply  perfect  correlation- 
neglecting  the  fact  that  ^  once  H2  desorption  begins). 

The  relative  yields  of  H^  a.:  !  HX  constitute  an  important  input  to  the  competitive  pairing 
model,  and  it  is  important  to  critically  examine  our  assumption  of  a  0.5  ML  saturation  coverage  for 
HX  and  to  assess  the  sensitivity  of  the  inferrt  ■'  oreferential  pairing  energies  to  this  assumption. 
Saturation  coverages  have  been  reported  for  a  number  of  hydrides  on  Ge(lOO)  and  Si(lOO),  and 
nearly  all  of  them  are  approximately  equal  to  0.5  ML.  Examples  include  H208^  and  H2S^^  on 
Ge(lOO)  and  H20,^*  H2S,*^  NH3,^0  C2H2,^*  and  C2H491  on  Si(lOO).  The  only  exception  to  the 
pattern  of  0.5  ML  saturation  coverages  of  small  hydrides  on  Si(lOO)  or  Ge(lOO)  of  which  we  are 
aware  of  is  HCl  on  Si(lOO),  whose  adsorption  has  been  reported  to  satuiate  at  0.25  ML.20  A 
saturation  coverage  of  0.5  ML  would  be  f'r.  ’ected  from  dissociative  adsorption  of  hydrides  AH^  as 
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H  and  AH^-i  on  surface  dimers  (Eq.  12),  or  from  di-a  bond  formation  in  the  case  of  C2H2  and 
C2H4.  Recent  careful  kinetic  uptake  measurements  on  some  of  these  systems  indicate  saturation 
coverages  of  slightly  less  than  0.5  ML,  viz.,  0.37  ML  for  C2H2  and  C2H4  on  Si(lOO).^*  Saturation 
coverages  slightly  less  than  one  molecule  per  surface  dimer  could  result  from  steric  repulsion,  from 
surface  defects,  or  from  occasional  dissociative  adsorption  events  across  dangling  bonds  on 
adjacent  dimers  rather  than  on  a  single  dimer.  The  latter  would  result  in  isolated  dangling  btmds 
that  may  be  relatively  inert  against  adsorption.^  On  Ge(  100),  saturation  coverages  of  both  Cl  and 
Br  have  been  reponed  as  1  ML,^'^  so  steric  repulsion  as  a  limiting  factor  in  HX  adsorption  seems 
unlikely.  The  use  of  relative  TPD  yields  to  determine  coverages,  as  in  this  study,  should  reduce  or 
eliminate  the  effect  of  non-adsorbing  defect  sites,  since  they  should  adsorb  neither  the  calibration 
molecule  (H2S  in  the  present  case)  nor  the  molecule  under  study  (HX).  Our  measured  yields  of 
0.30  and  0.31  ML  of  H2  desorbing  following  saturation  exposures  to  HCl  and  HBr,  respectively, 
excludes  a  saturation  coverage  of  0.25  ML.  Based  on  the  evidence  just  summarized,  it  appears 
most  likely  that  the  saturation  coverages  of  HCl  and  HBr  on  Ge(lOO)  are  slightly  less  than  0.5  ML. 

While  determination  of  the  precise  saturation  coverages  of  HCl  and  HBr  on  Ge(lOO)  will 
require  further  study,  we  find  that  our  conclusions  of  preferential  pairing  of  H+X  as  well  as  of  H+H 
are  not  sensitive  to  the  precise  saturation  coverage.  If  the  saturation  coverage  of  HCl  were  only  0.4 
ML,  the  observed  yield  of  0.30  ML  of  surface  hydrogen  desorbing  as  H2  would  imply  that  0.10 
ML  of  hydrogen  desorbs  as  HCl.  The  fraction /of  surface  hydrogen  desorbing  as  H2  would 
increase  slightly  from  0.75  at  saturation  coverage  to  0.79  at  an  initial  coverage  of  0.05  ML.  A  least- 
squares  fit  to  the  Tp  data  for  H2  and  HCl  (with  the  initial  coverages  scaled  downward  as 
appropriate)  using  the  modified  yields  produced  agreement  with  the  data  nearly  as  good  as  that  in 
Fig.  6,  with  f  in  the  range  of  0.75-0.74  for  coverages  between  0.04  and  0.4  ML.  The  parameters 
obtained  in  the  modified  fit  were  ehh  =  =  £Xx  =  2.85  kcal/mol,  f'a.HH  =  41.2  kcal/mol,  and 

^a,HX  =  38.4  kcal/mol.  The  inferred  pairing  energy  is  nearly  identical  to  that  obtained  assuming  0.5 
ML  saturation  coverage  (2.7  kcal/mol)  and  well  within  the  ±1  kcal/mol  uncertainty.  We  thus 
conclude  that  the  uncertainty  in  the  actual  saturation  coverage  of  HX  contributes  only  a  slight 
uncertainty  to  the  inferred  pairing  energies. 

Our  analysis  of  the  TPD  results  idealizes  the  surface  as  being  perfectly  flat.  However,  the 
sample  is  tilted  4-6*  from  (100)  and  consequently  7-10%  of  the  surface  atoms  are  at  step  edges. 
While  it  would  be  worthwhile  to  quantify  the  effects  of  steps  on  desorption  kinetics  from  (100)2x1 
semiconductor  surfaces,  we  argue  that  these  effects  are  minor  and  do  not  substantially  affect  our 
conclusions.  First,  in  the  only  directly  analogous  system  for  which  data  is  available,  Boland  found 
that  the  spatial  distribution  of  (paired)  hydrogen  atoms  desorbing  from  Si(  100)2x1  is  uncorrelated 
with  steps  or  defects.^’^^  Sites  vacated  by  desorbing  H2  following  brief  annealing  cycles  of 
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Si(l(X))2xl:H  were  imaged  by  STM.  At  high  hydrogen  coverages  desorption  events  occurred  at 
random  surface  positions,  while  at  lower  coverages  the  remaining  hydrogen  tended  to  form  one¬ 
dimensional  chains  due  to  effective  interactions  between  doubly-occupied  dimers.^*69.70,92  all 
coveragt  /:  examined,  however,  there  was  no  significant  correlation  between  desorption  sites  and 
step  edges.^’^  Second,  comparative  studies  of  adsorption  of  several  molecules,  including 
H20,84.89.93  H2S,89  CiUa,^  €2^2,^^  and  HCOOH,96  on  flat  and  stepped  Si(  100)2x1 
found  no  significant  effect  of  the  stepa  upon  adsorption.  By  detailed  balance,  desorption  cannot  take 
place  preferentially  at  steps  unless  adsorption  leads  to  preferential  population  of  step  sites.  Finally, 
if  H  and/or  X  adsorbed  preferentially  at  steps  these  sites  would  remain  saturated  at  coverages  above 
7-10%  (the  step  density)  and  could  not  affect  either  the  TPD  peak  temperature  nor  the  H2  versus 
HX  yield  at  higher  initial  coverages.  Since  the  data  used  to  infer  the  effective  pairing  energies  in  the 
model  were  for  coverages  above  0.1  monolayer  (cf.  Figs.  6, 7),  we  conclude  that  steps  make  at 
most  a  minor  contribution  to  the  observed  behavior. 

Having  concluded  that  preferential  pairing  of  H+X  occurs  just  as  for  H+H  on  Ge(lOO),  that  is, 
that  £hH  and  £Hx  (and,  probably,  exx)  are  greater  than  0,  we  now  turn  to  the  implications  of  the 
fitted  values  of  e.  If  the  n  bond  on  surface  Ge-Ge  dimers  provides  the  driving  force  for  preferential 
pairing,  then  one  might  expect  that  £hh.  Qtx.  and  exx  would  be  very  nearly  equal  to  chh  for 
hydrogen  adsorbed  alone,  or  approximately  5  kcal/mol.'^^  Qk  and  exx  might  be  reduced 
somewhat  because  of  steric  and/or  electrostatic  repulsion  between  the  paired  adsorbate  atoms, 
particularly  Cxx*  but  £^h  would  be  expected  to  be  the  same  as  when  no  halogen  atoms  were  present. 
Just  as  we  arc  able  to  conclude  that  the  pairing  energies  (2.75  ±  1  kcal/mol)  are  greater  than  zero 
(Figs.  8, 9),  it  is  also  clear  that  they  are  less  than  5  k  ;al/mol  (Fig.  8).  Several  explanations  are 
possible  why  the  effective  pairing  energies  are  less  than  5  kcal/mol,  which  cannot  be  resolved  at 
present.  One  possibility  is  that  £hh  is  close  to  5  kcal/mol  and  Chx  and  exx  arc  less  due  to  steric  and 
electrostatic  effects  and  that  the  fits  are  not  able  to  distinguish  this  case  from  the  equal-e  one. 
Experimentation  with  fits  with  various  pairing  energies  suggests  this  is  not  the  case,  however 
inferior  fits  are  obtained  with  einn  =  5  kcal/mol.  A  second  possibility  is  that  the  presence  of 
coadsorbed  halogen  atoms  does  in  fact  reduce  the  driving  force  for  pairing  (recall  that  excess 
halogen  atoms  remain  on  the  surface  during  a  TPD  ramp  while  H2  and  HX  desorb).  Such  a 
reduction  might  occur  because  the  more  electronegative  halogen  atoms  withdraw  electron  density 
from  the  surface,  which  decreases  the  population  of  the  7i  band  state.  Another  possibility  is  that  the 
adsorbed  halogen  atoms  cause  a  breakdown  of  quasiequilibium.  If  the  halogen  atoms  are  less 
mobile  than  adsorbed  H  atoms,  they  may  not  be  able  to  reach  the  equilibrium  extent  of  pairing 
during  a  TPD  ramp.  In  addition,  since  diffusion  along  the  dimer  rows  constitutes  a  quasi-one- 
dimensional  system,  they  may  also  block  diffusion  H  atoms  and  impede  pairing  of  H+H  and  H+X. 
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Further  study  will  be  necessary  to  determine  the  (Higin  of  the  quantitative  differences  between 
pairing  of  H+H  adsorbed  alone  and  in  the  HX  system. 

Estimates  of  the  heats  of  adsorption  of  HCl  and  HBr  can  be  obtained  from  both  the 
desorption  kinetics  and  from  molecular  bond  energies.  By  detailed  balance,  the  activation  energy 
for  desorption  of  HCl  or  HBr  should  be  equal  to  the  heat  of  adsorption  plus  the  activation  energy 
for  adsOTption.  The  near-unit  values  of  So  imply  that  the  latter  is  very  small,  as  discussed  above,  so 
the  heat  of  adsorption  is  approximately  40  kcal/mol  for  both  HCl  and  HBr  on  Ge(lOO). 

Independendent  estimates  of  the  heats  of  adsorption  based  on  molecular  bond  energies  are 
instructive  but  are  limited  by  the  paucity  of  accurate  thermochemical  data  for  germanium 
compounds.  Dissociative  adsorption  of  HX  yields  new  Ge-H  and  Ge-X  bonds  and  breaks  the  H-X 
bond  and  the  Ge-Ge  it  bond  on  the  dimer.  The  Ge-H  bond  energy  in  GeH4  has  been  measured  as 
82  kcal/mol.^^  However,  by  analogy  to  silicon,  replacement  of  the  other  Ge-H  bonds  in  GeH*  by 
bonds  to  other  germanium  atoms  on  Ge(l(X))  is  likely  to  reduce  the  Ge-H  bond  strength  on  the 
surface.  The  Si-H  bond  energy  in  ((CH3)3Si)3Si-H  (79  kcal/mol^^)  is  12%  weaker  than  that  in 
SiHa  (90  kcal/mol^^).  Assuming  that  a  similar  reduction  occurs  for  germanium,  we  arrive  at  an 
estimate  of  72  kcal/mol  for  the  Ge-H  bond  energy  on  a  (doubly-occupied)  Ge(lOO)  dimer.  Reliable 
germanium-halogen  bond  strengths  are  not  available,  and  for  the  surface  Ge-Cl  and  Ge-Br  bonds 
we  use  tabulated  average  bond  energies  of  81  and  66  kcal/mol,  respectively.*^  For  the  it  bond 
strength  in  Ge-Ge  surface  dimers  we  use  5  kcal/mol,  as  inferred  from  the  departure  of  the 
desorption  kinetics  from  first  order  at  low  coverage  for  hydrogen  adsorbed  alone.^  The  bond 
energies  of  H-Cl  and  H-Br  are  103  and  88  kcal/mol,  respectively.****  We  thus  arrive  at  72  +  81  - 
103  -  5  =  45  kcal/mol  for  the  heat  of  adsorption  of  HCl  on  Ge(lOO),  and  72  +  66  -  88  -  5  =  45 
kcal/mol  for  the  heat  of  adsorption  of  HBr. 

The  similarity  in  the  heat  of  adsorption  and  activation  energies  for  desorption  of  HQ  and  HBr 
on/from  Ge(lOO)  are  in  agreement  with  an  early  calorimetric  measurement  on  polycrystalline  Ge 
powder.  The  bond  strength  calculation  also  makes  clear  that  this  similarity  results  from  the 
weakening  of  the  surface  Ge-Br  bond  relative  to  that  of  Ge-Cl  being  very  nearly  equal  to  the 
difference  in  bond  energies  between  H-Cl  and  H-Br.  The  near  agreement  of  the  bond  strength 
estimates  with  the  activation  energies  for  desorption  suggests  that  the  bond  strength  estimates  are 
individually  reasonably  accurate.  More  accurate  determinations  will  be  necessary,  however,  for 
specific  inferences  to  be  made  about  individual  surface  bond  strengths. 

Preferential  pairing,  which  the  available  evidence  suggests  is  a  general  phenomenon  in 
adsorption  on  Si(lOO)  and  Ge(lOO),  will  affect  the  spatial  distribution  of  adsorbed  species  during 
film  growth  by  chemical  vapor  deposition  or  atomic  layer  epitaxy.  If  the  gas  phase  growth  species 
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is  a  stable  molecule,  two  adjacent  sites  on  a  surface  dimer  will  be  necessary  for  dissociative 
chemisorption  to  occur.  If  most  of  the  dangling  bonds  (and  adsoit>ed  species)  are  paired,  larger 
saturation  coverages  will  be  possible  with  shorter  exposures  than  would  be  the  case  is  the 
distribution  of  dangling  bonds  on  the  surface  were  random.  Near- unit  saturation  coverages  would 
be  particularly  helpful  for  atomic  layer  epitaxy  processes.  However,  the  3-5  kcal/nx)l  pairing 
energies  inferred  from  the  desorption  kinetics  are  small  enough  so  that  at  typical  growth 
temperatures  pairing  occurs  to  a  substantial  but  incomplete  extent.  Independent  measures  of 
preferential  pairing  could  provide  further  insight  into  this  phenomenon  and  its  usefulness  for  self- 
limiting  processing  chemistry  for  semiconductor  surfaces. 
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Fig.  1 .  TPD  spectnim  of  desorption  products  from  GeflOO)  following  a  saturation  dose  of  HCl. 
For  GeCl2  the  GeCl+  cracking  fraction  was  actually  detected. 


2.  TPD  spectrum  of  desorption  products  from  Ge(  100)  following  a  saturation  dose  of  HBr. 
For  GeBr2  the  GeBr^  cracking  fraction  was  actually  detected. 


m/e  =  36  amu  signal  (arb.  units)  m/e  =  2  amu  signal  (arb.  units) 


Fig.  4.  TPD  spectra  (2  K  s'*  heating  rate)  of  H2  and  HCl  desorption  products  from  Ge(lOO) 
following  various  doses  of  HCl.  The  initial  HCl  coverage  is  indicated  for  each  trace. 
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Fig.  5.  TPD  spectra  (2  K  s**  heating  rate)  of  H2  and  HBr  desorption  products  from  Ge(lOO) 
following  various  doses  of  HBr.  The  initial  HBr  coverage  is  indicated  for  each  trace. 


Fig.  6.  Dependence  of  TPD  peak  temperatures  for  H2  and  HCl  on  initial  coverage  of  HCl. 

Curves  throujh  data:  prediction  of  competitive  pairing  model  with  fHH  =  ^CI  =  fiCICI 
2.7  kcal/mol,  and  f'a.HH.  ^’a.HCl  =  41.4  and  38.1  kcal/mol,  respectively. 


Fig.  7.  Dependence  of  TPD  peak  temperatures  for  H2  and  HBr  on  initial  coverage  of  HBr. 
Curves  through  data:  prediction  of  competitive  pairing  model  with  £HH  =  %Br  = 
2.8  kcal/mol,  and  £a,HH.  ^aJlBr  =  41.7  and  38.5  kcal/mol,  respectively. 
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Fig.  8.  Predictions  of  competitive  pairing  model  with  various  values  of  Chh  =  ^HCl  =  fiClCI= 

compared  with  experimental  data  for  HCl.  For  £  =  5,  2.7,  and  1  kcal/mol  faJlH  was  41.7, 
41.4,  and  41.0  kcal/mol,  and  £a.HCl  was  38.6, 38.1,  and  37.7  kcal/mol,  respectively.  The 
latter  values  were  adjusted  slighdy  in  each  case  to  best  reproduce  the  TPD  peak 
temperatures  and  relative  yields  of  H2  and  HCl. 


9.  Predictions  of  competitive  pairing  model  with  different  assumed  pairing  energies 

compared  with  experimental  data  for  HCl.  Dashed  curves:  QiH  =  3.5  kcal/mol,  ChO  = 
ftici=  0.  £^a.HH  =  41.9  kcal/mol,  £a,HCl  =  37.9  kcal/mol.  Solid  curves:  £hh  =  ^HCl  = 
2.1  kcal/mol,  ficiCP  ^a,HH  =  41.1  kcal/mol,  £3, HCl  =  38.3  kcal/mol.  The  range  of 
predicted  values  of/,  the  fraction  of  hydrogen  desorbing  as  H2,  over  the  initial  coverage 
range  of  0.05  -  0.5  ML  is  given  in  each  case. 
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